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Abstract 

In 'arXiv:0908.2463 we computed the four-loop correction to a function 
depending on the 't Hooft couphng(s) that appears in the magnon disper- 
sion relation of the spin chains derived from single trace operators in = 6 
superconformal Chern-Simons theories. In this paper we give detailed de- 
scriptions of this calculation and the computation of the four-loop wrapping 
corrections for a length four operator in the 20 of SU{4), the i?-symmetry 
group for these theories. Here, we give all relevant Feynman diagrams and 
loop integrals explicitly, and also demonstrate the cancellation of double 
poles in the logarithm of the renormalization constant. 
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1 Introduction and summary of results 



The ABJM model is an A/" = 6 superconformal Chern-Simons (CS) theory coupled to 
matter [1]. The CS theory has gauge group U{N) x U{N), where the levels of the first 
and second gauge group are respectively given by k and —k. The matter multiplets 
transform in bifundamental representations of the two groups and in the fundamental 
or antifundamental representations of S'f/(4), the i?-symmetry group. This model also 
has a parity symmetry, where the spatial inversion is accompanied by an exchange of 
the two gauge groups, with a corresponding interchange of gauge representations for the 
matter fields. 

The gravity dual for the ABJM model is M-theory on AdS4 x S^/Z^ where the 
compact space S^/Z^ is equivalent to CP^ with a U{1) fibration. In the limit of large k 
the circle of the U{1) shrinks to vanishing size and the M theory reduces to type IIA 
string theory on AdS4 x CP^. 

This AdS/CFT correspondence for the ABJM model at large k closely parallels the 
more familiar correspondence relating SU{N) M' = 4 Super Yang-Mills (SYM) to type 
II B string theory on AdSs x S^. In both cases, the gauge theories are superconformal, 
with underlying supergroups PSU{2,2\A) and OSP{6\4) for the AT = 4 SYM and the 
A/" = 6 CS theories respectively. In particular, the spectra of the gauge theories should 
match the spectra of their string duals. On the gauge theory side the spectrum is given 
by the anomalous dimensions of the gauge invariant composite operators built from the 
elementary fields of the theory. In the planar limit, where N ^ oo the respective 't 
Hooft coupling constant A = Qym-^ ^r A/" = 4 SYM theory and A = y for the ABJM 
model is kept fixed. The relevant gauge invariant operators in this limit involve a single 
trace of the elementary adjoint fields in the case of A/" = 4 SYM theory, or a trace of an 
alternating product of bifundamental matter fields in the ABJM case. 

In both theories these composite operators mix under renormalization, and the 
anomalous dimensions are given as the eigenvalues of the respective dilatation operator, 
which is built from the pole part of the corresponding matrix-valued renormalization 
constant. This mixing of the composite operators can be mapped to a long range in- 
tegrable Hamiltonian acting on a spin chain. This allows one to find the anomalous 
dimension eigenvalues by solving appropriate Bethe ansatze [2H1]. The spin chains have 
ground states that correspond to the chiral primary operators in their respective gauge 
theories. These are operators whose dimensions are protected by supersymmetry and 
thus have zero anomalous dimension. A convenient choice for a spin chain ground state 
in A/" = 4 is tr(Z^), where Z is one of the complex scalar fields. In the ABJM model a 
convenient choice is ti^iY^Y^)^), where Y'^ is one of the scalars making up the 4 of the 
i?-symmetry group and is one of the scalars in the 4 [5H7]. This choice of a ground 
state breaks the supergroups to the subgroups SU{2\2) x SU{2\2) for A/" = 4 SYM and 
SU {2\2) for the ABJM model. The other single trace operators are then constructed 
by introducing magnons that change the fields in the chain. The magnons themselves 
transform in short representations of the respective unbroken subgroups. 

The presence of the SU {2\2) structures imposes severe constraints on the magnon 
dispersion relations of the respective Bethe ansatze. As was shown in [8J, these dispersion 



1 



relations must have the form 

E{p) = ^Jq^ + Ah^X) sin^f-Q, (1.1) 

where p is the momentum of the magnon on the spin chain, h?{X) is a function depending 
on the 't Hooft couphng, and A. Q is the i?-charge of the magnon. For A/" = 4 SYM the 
charge for a fundamental magnon is Q = 1, while in ABJM theory the charge is Q = \. 

The 't Hooft coupling enters the asymptotic Bethe equations only through this same 
function h^{\) [21110] • To the best of our knowledge, integrability makes no prediction 
for /i^(A), but there exist some alternative analyses not based on integrability [TT|[T2]. 
At weak coupling, it can be obtained order by order from perturbative calculations 
in the gauge theory, and the strong coupling expansion follows from the analysis of 
classical string states and their quantum corrections in the dual string theory. In A/" = 4 
SYM, all known results are consistent with /i^(A) = A/(47r^). However, in the ABJM 
model it is known at the level of two-loop perturbation theory that /i^(A) = A^ + (9(A^) 
[SlinilTS], while at large coupling we know by taking the BMN limit of type HA string 
theory on AdS4 x CP^ [6l[T3] and from one-loop string corrections [15] that /t^(A) = 
\\-^^+0{l) . In fact, since the perturbative expansion is even in A, the function 
/i^(A) should have a square root branch cut along the negative real axis. 

The calculation we present here also considers the ABJ modification of the ABJM 
model [IS], where the gauge group is generalized to U{M) x U{N), but with the levels 
of the gauge group kept at k and —k. In this case there are now two 't Hooft couplings 
defined as 

, M N , , 

but the superconformal group is still 05*^(614). It was shown at the two-loop level in 
the scalar sector [17J and in the full 05*^(614) sector [18j that A^ is replaced by AA, but 
otherwise the dilatation operator is the same, and is therefore still integrablejj If in the 
planar limit the integrability is to persist to higher loop orders for general values of A 
and A, then the only modification that can occur is in the function h"^. In general it will 
be a two parameter function, /i^(A, a), where we define A and a as 

A=v^, a=^^. (1.3) 

A 

Because the numbers of colour for the two gauge groups are different, the ABJ theory 
is not invariant under parity. Since a parity transformation exchanges the fundamental 
and antifundamental representations of the gauge groups, it has the effect of reversing 
the order of the fields inside the trace. In other words, it acts as a parity inversion on the 
spin chain. This spin-chain parity switches the odd and the even sites, or equivalently 
exchanges A with A. Clearly, the two loop result is invariant under paritjH, but at higher 
loops we should allow for different dispersion relations (11. ip for magnons at odd and 
even spin chain sites, with each having the form 



EoMip) = VQ2 + 4/i2(A,a)sin2|-Q , Ee.en(p) = ^odd(p)L . (1-4) 



^In the ABJM case, hints for a breakdown of intcgrabihty beyond the planar hmit have exphcitly 
been found in [IST. 

breakdown of parity invariance beyond the planar limit has been observed in |20) . 
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The ABJM case is captured by the special value a = 0, where the dispersion relations 
for odd and even sites automatically coincide. 

Up to four loops, the expansion of h'^{X,a) is then assumed to have the following 
form 

h\X,a) = \^ + X^hi{a). (1.5) 

The two loop contribution is independent of cr and the dispersion relation is obviously 
the same for even and odd magnons. At four loops, terms which are odd under a — — cr 
break parity and cause differences in the two dispersion relations. However, there could 
also be a dependent terms which preserve the parity. 

In this paper, we present in detail our computation of /i^(A,cr) to four-loop order. 
We already presented the final result in [21] 

In section [2] we show how it can be extracted from the four-loop dilatation operator 
in the flavour SU{2) x SU (2) subsector. The dilatation operator itself is reconstructed 
from an explicit perturbative four-loop renormalization of the respective composite op- 
erators as is summarized in section [3l The calculation can be simplified by considering 
only those graphs which give rise to non-trivial permutations of the flavour degrees 
of freedom. Nevertheless, we have to find and then evaluate explicitly more than one 
hundred diagrams plus their reflections where necessary. In section H] we present the 
classification of these diagrams, and their explicit evaluations. The dilatation operator 
and the result for /i^(A, a) is then extracted in section O 

In section E] we present our computation of the wrapping correction for length four 
operators in the SU{2) x SU{2) flavour sector. Our result matches prediction of Gromov, 
Kazakov and Vieira [23] which they obtained by means of the proposed F-system [231426] 
for the ABJM model. Finally, in section [7] we present our conclusions. 

Several details about the calculation have been delegated to various appendices. In 
appendix El we list the underlying component Feynman rules, as well as present use- 
ful effective Feynman rules which allow us to resolve the complicated spacetime-tensor 
structure of the graphs involving gauge bosons or fermions. In appendix [B] we sum- 
marize the flavour permutations which appear in the various Feynman graphs. Some 
properties of the appearing permutation structures are briefly discussed in appendix 
O The transformation properties of the Feynman graphs are discussed in appendix [Dl 
We then summarize in appendix [E] our two-loop computation of the scalar self-energy, 
which besides the known pole parts also yields the required finite parts. In appendix 
[F] an important sign is fixed by the two-loop renormalization of the permutation part 
of the six-scalar vertex. In appendix [Gl as a consistency check for our four-loop com- 
putation, we explicitly demonstrate the cancellation of double poles in the logarithm of 
the renormalization constant. The G-functions and triangle relations for an evaluation 
of the underlying loop integrals are summarized respectively in appendices [H] and [H 
Finally, appendix [J] contains the lists of the underlying several hundred two- three- and 
four-loop integrals which we need for the computation and for cross checks. 

■^Previous versions of this paper had an incorrect result for h^{X, a) due to several sign errors arising 
from inconsistent Feynman rules. A reexamination of the result was triggered by inconsistencies with 
an ongoing computation of h?{X,a) using the Af ~ 2 superspace formalism [35]. 
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2 Extraction of /i^(A, a) 



To compute the function h'^{X,a) in a field theory calculation, it is only necessary to 
consider a single scalar magnon either placed on the odd or the even sites of the spin 
chain. Thus, we may restrict our study to the SU{2) x SU{2) sector of the full OS'p(6|4) 
group, where the magnons in the first SU{2) live on the odd sites and those in the 
other SU{2) live on the even sites. The two different types of magnons do not interact 
with each other until the six-loop level pOj so for our purposes we may treat them as 
noninteracting. 

The dispersion relation for magnons at odd sites (11. 4p yields the following expansion 
in a power series in A 

Kddip) = 2A2(1 - cosp) + 2X\h4{a) - 3 + (4 - h^ia)) cosp - cos2p) + 0{X^) , 

(2.1) 

where we have set Q = 1/2 and we have reexpressed powers of sin | in terms of cosine. 
The exponentials which are contained in the cosine functions are powers of e*^, which 
is the eigenvalue of the shift operator that moves the magnon over by two sites on the 
spin chain. Hence, at the four-loop level there is a maximal shifting of four sites, which 
does not depend on h^. However, for shifts of two sites, there is an /i4 dependence. 

Without interactions between odd and even site magnons, the energy Eodd{p) is given 
as the eigenvalue of the dilatation operator, when acting on the (open) momentum / 
shift operator eigenstate 

L 

ijp = J2 e'''''iY^Yl)''Y^Yl(Y^Yl)^-^-^ (2.2) 

fe=0 

of a single magnon at odd sites. The momentum dependence in the expansion (12. ip 
is a result of acting with certain products of permutations on the above state, where 
a permutation exchanges the fields between two nearest neighbour either odd or even 
sites. Similar to the A/" = 4 SYM case, we introduce the permutation structures 

L 

{ai, 0-2, • • • , O-m} — P2j+ai 2j+ai+2 P2i+a2 2«+a2+2 • • • P2i+a„ 2j+a„+2 , (2.3) 

1=1 

where we identify L -|- i ~ i when we act on a cyclic state of length L. Some details 
about these structures and the permutation basis at four loops can be found in appendix 
[Cl To four loops, the dilatation operator expands as 

D = L + X^D2 + X^D^ia) + C(A^) . (2.4) 

At each order, the dilatation operator components decompose into a direct sum as 

Dk = -C^fc,odd "I" Dk,cven ~^ -Ofc mixed ; (2-5) 

where the individual parts act non-trivially respectively on odd and even sites only or 
mix odd and even sites. The permutation structures (12. 3 p which appear in the three 
parts therefore either carry only odd or even argument or both types of arguments. The 
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fact that even and odd site magnons are non-interacting to four-loop order immediately 
requires the absence of a mixing term Z)4,mixed involving the structures {1,2} and {2, 3} 
of the basis (]C.2p . We make an ansatz for i54,odd and apply it to the state (12.21) . The 



permutation structures thereby produce the following phase factors 

{}^1, {l}^e'P + e-'P , {1,3} -^e''^'+2e-'^ , {3, 1} e-''^' +2e'P . (2.6) 

This yields the eigenvalue as a function of the momentum p. Its comparison with the 
expansion of the single-magnon energy ( 12. ip then fixes the dilatation operator to 



D2,e.en = {} " {1} , 

D2Md = {}- {2} , 
A,oddM = (Kia) - 4){} + (6 - h,{a)){l} - {1, 3} - {3, 1} , 
I^4,even(^) = (/^4(-^) " 4){} + (6 - h^i-a)){2} - {2, 4} - {4, 2} . 



(2.7) 



We have added the part acting on even sites with the coefficients transformed by cr — t- —a 
as required by ( 11. 4p . The above expression yields zero energy for the ground state, since 
the sum of the coefficients of the different permutation structures sum up to zero. The 
above result allows us to find the function /i4(ct) by a direct perturbative computation of 
the dilatation operator. It is thereby sufficient to focus on D4^odd, and to compute only 
the contributions which involve non-trivial permutations. The coefficient of the identity 
part follows from the condition, that the ground state energy is zero. Computing the 
coefficients of the maximum shuffiing term and the absence of a part which mixes odd 
and even sites serve as some checks. These terms have also been computed in jTT] and 
the maximum shuffiing terms furthermore to six loops in [28|[29]. 



3 The dilatation operator from Feynman graphs 

The dilatation operator of the theory can be obtained from a perturbative computation 
of the one-point functions of the composite operators Oa- The appearing divergences 
due to quantum corrections require a renormalization of the composite operators as 

Ca.ren = 2j'Ob,ha.ve , Z = 1 + Z2 + \^Za, + . . . , (3.1) 

where the matrix Z cancels the appearing divergences and in general leads to mixing 
between the different composite operators. It is given by the negative of the sum of 
the divergences of the underlying Feynman graphs. In the formalism of dimensional 
reduction the spacetime dimension is reduced to 

D = 'i-2e. (3.2) 

The divergences then appear as inverse powers of e. The coupling constant A is accom- 
panied by the 't Hooft mass [i in the combination A/i^"^ to render the mass dimension of 
the loop integrals constant. The dilatation operator is then extracted as 

I? = /i-^lnZ(A/x2^£) . (3.3) 
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The logarithm thereby guarantees that all higher order poles in e cancel out, such that 
\yiZ only contains simple - poles. Inserting (13. ip . the expansion reads 



As a check for our calculation, in appendix [G] we demonstrate explicitly, that the double 
poles in the above four- loop contribution cancel out. The dilatation operator can thus 
also be obtained as 



Effectively, the above definition extracts the coefficient of the j pole at a given loop 
order L and multiplies it by a factor 2L. 

4 Four-loop diagrams 

In this section we first classify and then compute all logarithmically divergent simply 
connected four-loop Feynman diagrams with external scalar legs, which contribute to 
genuine flavour permutations. We neglect all contributions to the identity and trace 
operator in flavour space. These omissions are indicated by an arrow. We also stress 
that equalities between four-loop results hold up to irrelevant flnite contributions. 

We apply dimensional reduction, i.e. the algebra of the spacetime tensors is per- 
formed in exactly three dimensions. After that, the integrand of a given Feynman dia- 
gram only contains scalar products of momenta, and the integral can be dimensionally 
regularized in Z) = 3 — 2e. 

It is known [30] that this procedure is gauge invariant and hence preserves the 
Slavnov- Taylor identities at least to two loops. Thereby, in pure CS theory the statement 
holds for the pole part as well as for the flnite contributions, while in matter coupled 
CS theory only the pole part was considered. Instead in dimensional regularization, 
where the tensor algebra is also continued to D = 2 — 3£:, it is already observed for the 
divergent parts at two-loops in pure CS theory, that the Slavnov- Taylor identities are 
not fulfllled. It has not explicitly been checked whether dimensional reduction in CS 
theories is gauge invariant to four-loop order, but a breakdown appears very unlikely. 

4.1 Classification of Feynman diagrams 

We have to flnd all simply-connected subdiagrams which, when attached to the compos- 
ite operator, contribute to the renormalization at four-loops. Such diagrams at L loops 
with Vi elementary vertices with i legs fulflU 




(3.4) 



V = \im 2e\-^\nZ{\e) 



e^O dA 



(3.5) 




(4.1) 



At L = 4 loops, we obtain the following combinations of vertices 



type multiplicity 



Vis 2 1 1 1 
021043210 
002402468 



(4.2) 
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We only want to compute the diagrams with scalar external lines. They should also con- 
tribute to genuine flavour permutations. In this case the only relevant classes (Ve, V4, V3) 
are given by (2,0,0), (1,2,0), (1,1,2), (1,0,4) (0,4,0), (0,3,2). The members of each 
class in terms of the combinations of the vertices (lA.Sp of the theory are given in table 
[TJ All other combinations either only contribute to the identity or subtraces in flavour 



(^6, ^4, V2) 




V(yty)3 


1 fp'I tp ' 


11' yjyj ' 

V J. 

*^yty^t^ 


^AYAYt 


V A 4 Wt 

^AAYty 


^AipH 


^AY'iY 


Va3 


f2 OU 
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(1,2,0)| 
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1 





2 
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1 
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a 


1 





1 








1 


1 





(1,1,2) 1 


b 


1 








1 








2 





c 


1 











1 





2 







d 


1 











1 





1 


1 




a 


1 














2 


2 





(1,0,4) 1 


b 


1 

















4 





c 


1 

















3 


1 




d 


1 

















2 


2 


(0,4,0)| 


a 





2 


2 

















b 








4 



















a 








3 








2 








(0,3,2)| 


b 








3 








1 


1 







c 








3 











2 






Table 1: Combination of vertices which lead to planar four- loop diagrams involving 
permutations. The gauge-ghost vertices are not presented, since it is clear that they 
contribute whenever an gauge boson loop with cubic gauge vertices is present. 

space, or they contain tadpoles. 

The diagrams which can be built up from the combinations in tabled] will be ordered 
according to their flavour structure, i.e. we will classify them with respect to their 
number of six-scalar vertices and of quartic vertices involving fermions. 

4.2 Diagrams involving two six- scalar vertices 

In the following we evaluate all Feynman diagrams with two six-scalar vertices, which 
contribute to the coefficient of genuine flavour permutations in the dilatation operator. 
They build up the class (2, 0, 0) in tabled! 

Using the flavour permutation structure flB.2l) of the six-scalar vertex, the diagrams 
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which lead to non-trivial flavour permutations are given by 



So 



5*6 




M'N'h{{l,3}-{l}) 
(-i + ^)({l,3}-{l}), 
M'N'h^i{l,2} + {2,3}-{l}-{2}) 



16 

(Any 



(-4^ + ^)({l,2} + {2,3}-{l}-{2}), 




^-^M3iV/4bbb7({l} + {2}) 



1 

'4^ 



16 86 



(4.3) 



({l} + {2}) 



M2iV2(-K(/2)')(-{l}) 



(AA)^ 1 

16 £2 



{1} 



^7 



— > 



(470 



16 V2£2 ^ 



The explicit expressions for the appearing integrals can be found in (jJ.SSp . As explained 
at the beginning of this section, we have neglected all contributions to the trace and 
identity operator. These truncations are indicated by the arrows. We have also not 
evaluated the remaining scalar diagrams 




So 



s 



10 



(4.4) 

which only contribute to the trace or identity operator. 

Along with the set of diagrams in (14. 3p . we also need the reflected diagrams of 5*2, S'4, 
5*7, which are easily obtained by using the reflection operator R in fID.Sp . The diagrams 
S'4 and 5*5 contain permutations between odd as well as even sites. To compute the 
dilatation operator for odd sites and to check the vanishing of the terms which mix odd 
and even sites, we therefore have to consider the respective odd and mixed contributions 
from the diagrams 8(5*4) and S(S'5), which are shifted with the shift operator S in (ID.2p 
by one site. We thereby keep only those permutation structures, which have an odd 
entry as their first entry, which is indicated by the suffix {1, . . . }. The sum of the 
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relevant terms then reads 



S = S2 + RiS^) + [S^ + RiS^) + S(54) + S(R(54)) + ^5 + HS5)i^,^^^^y 

+ ^6 + ^7 + R(^7) 



AA 
16 



AA -T^ + T ({1>3} + {3,1} + 2{1,2})+ -- - 12--7r2 {1} 



2^2 



+ (A-A)^k{l} 



(4.5) 



4.3 Diagrams involving a single six-scalar vertex 

In the following we evaluate all Feynman diagrams with a single six-scalar vertex and 
additional interactions, which promote the diagrams to four loops. Their flavour per- 
mutation structure is the same as for the simple two-loop diagram which contains the 
six-scalar vertex, if the additional interactions inside the diagram form flavour-neutral 
interactions, i.e. they do not alter the flavour flow. These diagrams build up the classes 
(1,2,0), (1,1,2) and (1,0,4) in table [H 



4.3.1 Diagrams involving flavour- neutral next-to-nearest neighbour inter- 
actions 

Flavour-neutral next-to-nearest neighbour interactions can only be realized with gauge 
boson exchange and at least two cubic gauge scalar vertices at the flrst and last leg. 
The relevant diagrams involve the following substructure 



(4.6) 



They have to contain a six-scalar vertex which has the flavour permutation structure 
given in (lB.2p . and hence contains a simple permutation. The diagram involving the 
first of the above structures then spans the entire subclass (1, 1, 2)^ in table [H The 
second and third one belong to (1, 0, 4)^ . 

Using the effective Feynman rules (1A.7I1 . the only non- vanishing contributions after 
a convenient choice of the external momenta are given by 



^vl 




(4vr)\,2..2,. , .r ^M._(AA)^l^o 



^^M^iV^(/,,, + 2/,,,){l} g 



-2--U1}. (4.7) 



The corresponding four- loop integrals are explicitly given in (IJ.SOp . 



4.3.2 Diagrams involving flavour-neutral nearest-neighbour interactions 

An interaction between three scalar field lines 









+ 




(4. 



appears in diagrams of the type 



(4.9) 



These diagrams belong to the subclasses (1,1, 2)c and (1,0, 4)c in tabled! They are 
all finite, and in particular can be made vanish, if the external momenta are chosen 
conveniently, carefully avoiding the appearance of IR divergences. 

The flavour-neutral nearest-neighbour interactions involve the following diagrams 



M 



+ 



(4.10) 

In the flrst term, which contains a fermion loop, the flavour flows from the upper to 
the respective lower external line as in all the other diagrams. Several of the above 
subdiagrams vanish identically, as follows from the effective Feynman rules flA.Qp . 

The above subdiagrams appear in four-loop diagrams which also contains a six-scalar 
vertex, and hence according to (IB. 21) contributes to a simple permutation. They belong 
to all subclasses of (1,2,0), (1,1,2), (1,0,4) in tabled! 

Evaluating all the diagrams in which the six-scalar vertex is connected by the struc- 
tures fl4.10p to one of its nearest neighbours, we obtain with the effective Feynman rules 



Snl 



S, 



n2 



16 

(470 
k^ 



2/42bbb2 — -^Ii2bh4ab + Ii22do7adbd ) {1} 



(4.11) 



^ -^M3iV^/42bbbl{l} = ^ 

2 io 



1 1 \ r . 



In S-nii the contribution from the fermion loop has to be considered, while in 3^2 the cor- 
responding integral does not have an overall divergence. Only those diagrams in fl4.10p 
with gauge bosons contribute, which do not have a cubic scalar gauge vertex at one of 
the lines which is an external line. This can be seen if we set the external momentum at 
the corresponding line to zero. It is an IR safe choice and makes the corresponding con- 
tributions vanish. The reflected copies are found by applying the reflection operator R in 
( ID.3P that preserves the action of the permutations on respectively odd and even sites. 
For all inserted substructures f l4.10p which lead to logarithmically divergent diagrams, 
one flnds that R acts by simply exchanging A 4-> A. The sum of these contributions is 
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then given by 



{I}- 



(4.12) 



'S'nn — 'S'nl + R(S'nl) + 3^2 + R(5'n2) 

The flavour-neutral interactions (I4.10p can also involve two external lines of the 
six-scalar vertex. In this case, the relevant contributions read 



v2 




1 (An) 



2 



-M^iV K(J22a)K(/2){1} 



16 4^2 



{I}- 



(4.13) 



Only the diagram with two quartic scalar gauge vertices in f l4.10p contributes in this 
case. The overall pole part of the four-loop integral is given as the negative of the pole 
parts of two two- loop integrals, which are given in ( ]J.2p and ( \J.3\i . 

The flavour-neutral interactions fl4.10p can also involve one external and one internal 
line of the six-scalar vertex. The relevant contributions read 



v3 



— )■ 



(470 



-MN'H - -/, 
37r2 



42bb4afe 



+ 1 



422ch7adbd 



-21 



42bbb2 



{1} 



(4.14) 



1 

16 V 8i2 8e 

In this case, from fl4.10p only the first diagram with a fermion loop and the ones with 
a quartic gauge-scalar vertex at the external line contribute. The relevant integrals are 
explicitly given in flXiOj) . (1X43]) . flX39|) 

Finally, the diagrams with the flavour-neutral interaction f l4.10p involving two inter- 
nal lines of the six-scalar vertex are evaluated as 



-^-^ {1} 



v4 




AA 
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K{G{2 - 3A, l){M^N\h^^ + 2h^^) 

+ M3iV(/3gb - K(/22a)/2 + 2/3gt 

+ /3g, + 2/3gc + /3tb))){l} 

1 1 7 



(4.15) 
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4e2 eV 12' 

This is the most complicated contribution in this class, since it involves the maximum 
number of different contributions, that individually require intricate evaluations. The 
above result is expressed in terms of three-loop integrals, which involve the nearest- 
neighbour interaction fl4.1Up and the two additional loops built by the four attached 
propagators. The integrals are explicitly given in f |j.37p . 

We still have to include the reflections of the above contributions, given by acting 
with R in (]D.3p . Again, this corresponds to simply exchanging A O A. The sum of the 
above contributions is given by 

S'v = + 5'v2 + R(»S'v2) + *S'v3 + R(S'v3) + >S'v4 + R(S'v4) 



AA 
16 



AA 



5 1 / 41 . 

— + - 20 r 

Ae^ e\ 12 
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1/^ 23 . 

- 1 TT^ 

e\ 24 



{I}- 
(4.16) 



4.3.3 Diagrams involving the scalar two-loop self-energy 

The two-loop self-energy correction of the scalar field also appears as a sum of subdia- 
grams. It reads 
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(4.17) 

To calculate the divergences of the corresponding four-loop diagrams, we need the di- 
vergent and finite parts of the two-loop self-energy. The individual contributions are 
evaluated in appendix [El Apart from the Wick rotation (the result has to be multiplied 
by a factor i"^ = —1), the amputated scalar self-energy contribution becomes 



4 \2e 
(A-A)V l 
4 \Ae 



7t' + 3[ 
2 V 3 



7 + In 47r 



+ 



Y + 2(3-7 + ln47r) 



(4.18) 



-l+2e 



where the last propagator factor on the r.h.s. captures the momentum dependence. Its 
weight label indicates the exponent of ^, where p is the external momentum. The pole 
part of the above result coincides with the result in [T7] . 

The following diagrams contain the two-loop self-energy and contribute to the renor- 



malization of the composite operators (factor 
and the six-scalar vertex are included) 



i from the four scalar propagators 



Ssl 




5s2 







^ -^^^MNK{J:Y)K{h){l} 

aa|, + (a-a)^^]{i} 
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" 16 
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16 



MNK{^Yh{l + 26) - K(Sy)/2){l} 

3 



ll + -7r^ 
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1 1 
+ - 
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{I}- 
(4.19) 
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They belong to all subclasses of (1, 2, 0), (1, 1, 2), (1, 0, 4) in table[TJ The combination of 
the above diagrams, which enters the renormalization of the composite operator, reads 



Ss 



-Ssi + 35*82 
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16 
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(4.20) 

The sum of all contributions which involve a six-scalar vertex and a flavour-neutral 
four-loop completion is given by 
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(4.21) 



4.4 Diagrams involving four scalar-fermion vertices 

In the following we evaluate all diagrams of the class (0, 4, 0) of table [H The fermions 
have to form a single loop, which in this case is a square, to contribute to flavour 
permutations. The only four-loop diagrams of this class in which four neighboured 
scalar lines interact are given by 



— > 



1 {4ny 



s2 



8 

16 £ 



■M^N^h2bbMl, 2} + {2, 3} - {1} - {2}) 
({1,2} + {2, 3} -{1}- {2}), 



(4.22) 
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--^^M'^Nh2hhide{id. and tr.) 
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16 32^' 



fid. and tr.l 



where the r.h.s. follow immediately, using the effective Feynman rule for the fermion 
square flA.16p . The resulting integrals are explicitly given in f|j.42p and f|j.40p . and the 
flavour permutation parts are taken from flB.3p . The second of the above diagrams only 
contributes to the identity and trace part of the dilatation operator and hence is not 
considered here. 

Besides the above diagrams, there are the following diagrams in which only three 



13 



neighbouring lines interact. They can contain a bubble and are given by 
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(4.23) 

where the r.h.s. follow immediately, using the effective Feynman rules for the fermion 
squares with bubble flA.lTp . The resulting integrals are explicitly given in f|j.40p . f|j.42p . 
flJ.411) . and the flavour permutation part is taken from fIB.Sp . The diagrams with four 
scalar-fermion vertices can also contain triangles and read 



-Pstl 
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1 (47r) 
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16 e 
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(4.24) 
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where the r.h.s. follow immediately, using the last of the effective Feynman rules in 
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(lA.17p . The resulting integrals are explicitly given in (]J.45|) and f|j.40|) . and the flavour 
permutation part is taken from (IB. 31) . 

Summing up the above contributions, considering also the contribution from the 
diagram Fgi shifted by one site with the shift operator S in (ID.2p and allowing for 
factors of two that come from reflections of the diagrams Fghi to Fghe and -Fst2) "we 
obtain 

-^s = [-^sl + S(Fsi)]{i_... } + 2(Fsbl + -Fsb2 + -^sbS + -^sb4 + Fsh5 + -^sbe) + -^stl + 2Fst2 



16 



4 

-{1,2} 



(4.25) 

where we have only kept those permutation structures, which have an odd entry as their 
first entry as indicated by the suffix { 1 , . . . } . 



4.5 Diagrams involving three scalar-fermion vertices 

In the following we evaluate all diagrams of the class (0, 3, 2) of table [T] 

The four-loop diagrams in which three scalar-fermion vertices form a fermion triangle 
also involve a single gauge boson propagator. According to the flavour permutation 
structures in f lB.3P they contribute to the coefficient of a single permutation. The gauge 
boson propagator can end on two cubic gauge-fermion vertices. The corresponding 
diagrams build up the class (0, 3, 2)^ and are given by 
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(4.26) 
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'42bb3ad 



){1} 



A=^A 1 



where in the first and third line the equalities between the graphs are to be understood 
for the pole parts only. The r.h.s. of the diagrams involving the fermionic self-energy 
follow immediately from the effective Feynman rule f lA.12|) . The r.h.s. of the other 
diagrams are obtained from the effective Feynman rule (]A.13p . The resulting integrals 
are exphcitly given in (lJ.4ip . f |j.40p . ( ]J.43p . 
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The gauge boson propagator can also end on one cubic gauge fermion vertex and 
one cubic gauge-scalar vertex. The corresponding diagrams build up the class (0, 3, 2)^ 
and are given by@ 



tv3 




1 (AttY 

^ -^4^5^^'^'(^422qAaM " /422qAf,ad) (-8{1}) 
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16 3 ^ ' 
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(4.27) 
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TT 
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The r.h.s. of the diagrams are obtained from the effective Feynman rules (1A.15I1 . The 
resulting integrals are explicitly given in (IJ.41I) . (]J.40I) . (1J.43I) . ( 1J.45I1 and (1J.46I) . 

The remaining diagrams, in which the gauge boson propagator is attached to the 
diagram via two cubic vertices involving scalars build up the class (0, 3, 2)c and can be 
made vanish by a convenient choice of the external momenta. 

The above results depend on a sign z = ±1 which is not uniquely given in the 
literature (TJEl] . In appendix |F] we fix z by computing the renormalization of the six- 
scalar vertex at two loops. It should vanish for the correct sign choice of z to ensure 
superconformal invariance. This yields z = 1. 

We also have to consider the reflected diagrams. Their contributions are obtained 
by applying the reflection operator R in (lD.3p . We obtain for the sum of the respective 
diagrams 



Ft = 2(Ft,i + R(Ft,i)) + Ft,2 + R(Fts2) + 2(Ftvi + R(Ftvi)) 

+ Ftv2 + R(-^tv2) + 2(Ftv3 + -^tv4 + Ftv^) 



-)■ z- 



XX 
16 



XX{-^ + - 



(4.28) 



12 



7r 



TT 



+ (A-A)^- {1} 



The presence of the reflected diagrams which do not differ by a sign from the original 
diagrams, but only by an exchange A ^ A thereby guarantees that the result only 
depends quadratically on the difference of the couplings and hence on the parameter a 
defined in (11. 3p . 



^Because of inconsistent Feynman rules these three diagrams had a sign error in previous versions 
of this paper. 
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5 Result 



Fixing 2; = 1 as determined in appendix [Fl the part of Z involving non-trivial permu- 
tations is given as the negative of the sums of (14. 5p . (I4.2ip . (I4.25p . and (14.280 . 

-^)({1,3} + {3,1}) + 1{1,2} (5.1) 

By the arrow we indicate that we have neglected all contributions to the identity and 
trace operator in flavour space, and also that we have only written half of all contribu- 
tions, namely those with permutation structures with an odd entry as their first entry. In 
the above result the simple poles of the structure {1,2}, which couples magnons at odd 
and even sites, between the scalar contribution (14. 5 p and the one with fermion square 
(I4.25P cancel out. Only a double pole part is left. The neglected terms which contain 
permutations that have an even entry as their first entry can easily be reconstructed by 
acting with the shift operator S in (ID.2p on the above result. 

In appendix O we explicitly demonstrate that all the double poles cancel against 
each other when we add the two-loop contribution and take the logarithm as in (13. 4p . 

According to the prescription (13. 5p . the dilatation operator for odd sites is the co- 
efficient of the ^ pole terms in (15. ip multiplied by 8. We must still add the neglected 
identity part, which we fix by demanding that the ground state has zero eigenvalue. We 
then obtain 

I^4,odd = -(4 + 4C(2) + a\{2)){} + (6 + 4C(2) + o\{2)){\} - {1, 3} - {3, 1} , 

(5.2) 

where C(2) = \. By comparing this result with (12. 7p . we immediately see that the 
coefficients of the maximal shuffiing terms match. The four-loop term of the function 
/i^(A, cr) in (II. 5p is read off as 

/i4(a) = -(4 + a2)C(2). (5.3) 

6 Four-loop wrapping interactions 

To find the correct anomalous dimensions of a length four state at four loops, we have 
to consider the wrapping interactions [32l|33]. Recently, Gromov, Kazakov and Vieira 
have made a prediction for the four- loop wrapping contribution in terms of /i^(A) for 
scalar operators in the 20 representation of S'f/(4) [23] • The highest weight of this 
representation is in the SU{2) x SU{2) sector. These authors find their result by applying 
the thermodynamic Bethe ansatz (TBA), first introduced for the original AdS/CFT 
correspondence in |3lll35], using the predicted asymptotic Bethe equations [10]. In 
particular, the TBA is formulated in terms of a F-system [23H2S], a series of difference 
equations, that lend themselves to an efficient order by order solution. 
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Following the strategy of [36l[37], we first subtract from the asymptotic dilatation 
operator (12 .Zp the range five interactions (given by 5*2 and the reflected diagram). This 
leaves for D4^^odd in (15. 2 p 



D 



sub 
4,odd 



(4 + 4C(2) + a^C(2))({l}-{}). 



(6.1) 



We then have to add the following wrapping diagrams which contribute to genuine 
flavour permutations 
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We then sum the above diagrams as 



W = Wi + W2 + 2W3 + Wi + W5 



(AA) 
16 



2 r 



2e- 



TT 



{1} 



(6.3) 



where we include a factor of two for W3 since its reflection cannot be mapped to itself by 
a cyclic rotation. To this sum we add the identity terms needed to give zero wrapping 
for a chiral primary. Taking the negative of the sum W, extracting the residues of - 
and multiplying by 8 as required by (13. Sp . we find that the wrapping contribution to 
the odd site dilatation operator is given by 



D 



4,odd 



-(4-2C(2))({l}-{}). 



(6.4) 
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The 20 in the SU{2) x SU{2) sector is a completely antisymmetric state, so the 
eigenvalue of {1} — {} is —4, and we have to double this to take into account also 
the contribution from the even site dilatation operator -D^even- Hence the four-loop 
wrapping contribution to the anomalous dimension is 

7^0 = (32 - 16C(2))A^ , (6.5) 

which agrees with the GKV Y-system prediction [23] • It does not depend on o"^, which 
is expected from the topology of the diagrams and the fact that the two-loop result does 
not depend on a^. Adding to DJ, we find that the anomalous dimension of the 20 
is given by 

720 = 4 + 8P - 8(6 + ct2)^(2)A\ (6.6) 

where we have included the lower orders [5] . Intriguingly, the rational terms of (16. ip and 
(16.41) cancel, such that the four-loop contribution only consists of terms with maximum 
transcendentality. 

7 Conclusions 

Our four-loop calculation of the dilatation operator shows that the unknown function 
/i^(A, a) with the expansion (II. 5p has the form 

h'^{X,a) = X'^ + X^h4{a) , hi{a) = h4 + a^hi,a . (7.1) 

It is an even function of a, and hence to four loops there are no parity breaking effects 
which could lead to different dispersion relations for magnons at odd and even sites as 
in (II. 4p . The coefficients are explicitly given by 

/i4 = -4C(2) ^ -6.58 , /i4,. = -C(2) . (7.2) 

The negative sign for is sensible, as this will dampen the quadratic behaviour found 
at small A to the linear behaviour at large A. Interestingly, the rational contributions to 
/i4(o") cancel, leaving a maximally transcendental result. Inserting the definitions (II. 3p . 
it becomes obvious that due to the relative coefficient between and h^^a the four loop 
term depends on the coupling constants only via a factor (A + A)"^. 

The calculation is rather complicated and contains large numbers of contributions. 
We therefore have performed several cross-checks on the integral tables, and we have 
reproduced the known results at two loops. Furthermore, by checking the cancellation 
of double poles of the four-loop contributions in the logarithm of the renormalization 
constant, we have another check for our final result Jf] Finally, also the four- loop anoma- 
lous dimension of the non-protected length four state matches the result obtained from 
the F-system. 

^The same result is found using the Af = 2 superspace formalism [22] , 
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A Action and Feynman rules 

In three dimensional spacetime with metric rjf^i, — diag(— , +, +) The antisymmetric 
tensor is normalized as e^"^^ = 1. The product of 7-matrices having index positions 
{'yfi)^ is given by 

where z = ±1 represents a sign. The Clifford algebra is given by 



(A.l) 



1 
1 



The kinetic and interaction parts of the action are given by 



Svm — 



47r 



d^xtr 



'S'int = ^ I d^x tr 



-ze^^-'iAo^A^A.-A^A^A^) 

- iAf^Y^d^Yl - iAf^Yld^Y^ + 2YIA^Y^A'' 



(A.2) 



A-nAj^Y^Y A^AijY Y^ 



- iI^^^'Mb + A^ilj^'^^^'il^B - lA^ [c,d^c*] - %A^ %d^c*\ 

+ ^Y^yIy^Y^Y^yMsES^ + S^6E6E - Q6^^6E6E + AS^S^SE) 

- '-{YlY^i^^^i^n - i^ni^^^Y^Yi){SiSS - 2S^SE) 



+ i^e^'^'^'^YX^PBYl^Pn - ^eABcnY^^l^^^Y^ ^P^'' 



(A.3) 
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A.l Component Feynman rules 

We obtain the propagators for the gauge, scalar, fermion and ghost fields from S'kin ai 

{A^{p)Ap{-p)) = -(^(p)i^(-p)) = -ii-(e,^^p7 _ ,^P^ 



2p^ \ p 



5| 

5^ 



2 



(A.4) 



P 



{c{p)c*{-p)) = -i\ 



p2 



where diagonality in the gauge group indices and a factor for each propagator have 
been suppressed, and we will fix the gauge as C = in the following. 

The vertices are obtained by taking the functional derivatives of i times the action 



previous version of this paper had inconsistent Feynman rules. The propagators in (|A.4p and 
an overall sign in the scalar-gauge three-point vertices in (jA.Sp have been corrected. As a consequence 
this ahers the overaU signs in (|A.6p . (jA.lip . (|A.14p and (|A.15j) . 
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w.r.t. the corresponding fields. We obtain for tlie cubic vertices 
= Va3 = - = -2e"^^ tr (T" ) , 




/3 
y 



yt 



y 



c 



y 



yt 
y A* 



-'^ = K^tAv- = -^tr {B,T''B')Y5^ , 



= Vacc = Vac," = Pc* tr {T^ [T^T'=] ) , 



= VayAy^ = 2^tr {T'^B'BJ'y^r,^^'' , 



= I^Ayyt = T^iiyty = -^tr {T''B'B,T'')S^r'' , 



yt 

y V't 



yt 

y V'''' 



^VVtyyt = -VYiYi,H = -^tr {BaB''B,B''){S^S^ - 26^6^) , 

Vy^^Yf^ = -Vy.^YH = ^ ( {B„B,B,,B,)e^''^'' + 3 perm, ([^-^^j)) 
y ^ ,ABCz?(^j, {BaB,B,Ba) - tr {BaBaB.B,)) , 

Y Y"! Y 



A 

yt y yt 



F(yyt)3 



^ ( tr {B^B,B<^BaB^Bf) {F^gi + F^i^ + F™) + 11 perm. ) , 

(A.5) 

where the momenta carry as suffix the respective field line label at which they enter 
the vertex. We have suppressed a factor ^ for each vertex. The dependence on the 
coupling constant can be easily restored at the end of the calculation. 
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A. 2 Effective Feynman rules 



The gauge propagator in flA.4l) and the cubic gauge vertex in (lA.SP contain e tensors 
which appear in the numerators and are contracted with each other or with loop and 
external momenta. The same happens in presence of fermion fields. Dimensional re- 
duction requires, that the e-tensors are reduced to scalar products in strictly D = 3 
dimensions, before the integral is dimensionally regularized by switching to D = 3 — 25 
dimensions. It turns out to be advantageous to introduce effective Feynman rules, in 
which the corresponding tensors have already been reduced, such that all gauge bosons 
and fermions appear as scalar propagators with momenta in their numerators. 

In the following all substructures are given with amputated external legs. All pref- 
actors (apart from powers ^) of the propagators (1A.4I) and vertices (lA.Sp . symmetry 
factors (like —1 for a fermion loop) and factors for internal flavour loops are included 
in the prefactors. The corresponding factors for the propagators which have to be at- 
tached at the external legs are not included when the subdiagram is replaced by its 
scalar representative. 

Since the free vector indices of the gauge boson propagators are perpendicular to the 
momentum, we can deflne the following effective vertices (extracting the factors of i) 
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yt 
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■f^ -> — 2i 



yt 
yt 



~2i 



yt 
yt 



(A.6) 
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Reexpressing the appearing products of two e-tensors in the numerators of the loop 
integrals in terms of the metric, the flavour-neutral interactions via gauge bosons be- 
tween three scalar fleld lines are simplifled as 







= 2i| 
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(A.7) 



In the above relations factors —i for external propagators have not been included in the 
prefactor. The internal scalar propagators on the r.h.s. have no non-trivial prefactors. 

For flavour-neutral interactions between two scalar fleld lines involving the gauge 
bosons we need the cubic gauge vertex with propagators. It is reexpressed in terms of 
the metric as follows 



/3 - VafiV^/aVpft + VanV-yaVfip + Va/SV^/pVaK " Va-yVfipVan ■ 

(A.8) 
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The effective Fcynman rules then become 
1 



+ 4" 



0, 



— 



it /i 



(A.9) 

where we have omitted simple factors of either M ot N from the colour running in the 
loop. In the last diagram we have summed up the scalar, fermion, gauge boson and 
ghost contribution to the gauge boson self-energy. 

For the fermion loop wc find (considering also the factors 4 for the fiavour contraction, 
—1 for the fermion loop and another —1 for changing the direction of one momentum 
factor in the numerator) 

(A.IO) 



The fermion fields get finite self-energy corrections at one loop. We can summarize the 
one-loop fermionic self-energy as 



(A.11) 



where we have also included the non-trivial colour factor. The above substructure has to 
be considered in certain diagrams with a fermion loop. The relevant effective Feynman 
rule for the substructure is 



Y 



z-{M -N)M 




(A.12) 
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The constant z = ±1 thereby parameterizes a sign discrepancy in the hterature. For a 
gauge interaction between two sites of the triangle we obtain with — 1 the following 
rule 



Y f{ Y 

y 



—z- 






(A.13) 



The above rule simplifies, if the external momentum at one of the three scalar-fermion 
interactions can be set to zero. Two of the first three contributions on the r.h.s. then 
cancel against each other, and the last term vanishes. This is the case in all required 
diagrams at four-loops. 
The subdiagram 
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(A.14) 

is only needed in two special cases, in which either pi = or pa = p^. The effective 
Feynman rules respectively read 




(A.15) 



The substructures involving only scalar-fermion vertices lead to the following effective 
Feynman rule for a fermion square 
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(A.16) 
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We then immediately obtain the rules for the remaining relevant subdiagrams 



Y 




(A. 17) 



yt 



As in the case of a fermion triangle, also in this case the rules simplify if the external 
momentum entering one of the vertices which participates only in a single loop of the 
full four-loop diagram is set to zero. 



B Flavour structures 

Besides the momentum-loop and colour factors of the diagrams, we also have to work out 
their flavour structures. Therefore, in the following we focus on the flavour non-neutral 
vertices and their combinations and compute the flavour flow through the different 
structures. For our computation we only need the terms which involve non-trivial flavour 
permutations without any subtraces. 

The flavour structure of the six-scalar vertex explicitly reads 



BDF I T?DFB 
+ ^CEA 



^FBD\ 
EAC ) 



(B.l) 



It can be visualized as 



X 



( + ) 





2 ) ( +4 >< 



(B.2) 

The last term indicates the non-trivial permutation of the flavour at two next-to-nearest 
neighboured sites. It appears with prefactor one. The other terms can also lead to non- 
trivial permutations without subtraces when two six-scalar vertices are combined as in 

(iD. 

We have to work out the flavour structures of the following contractions of vertices 



26 



involving fermions 



I 



■'< 



+ 4 



+ 4 



+ 4 +4 
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-2)|( 



r 
r 
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■8 ) ( +8 >< ' (B-3) 
+ 20 ) ( -16 >|< , 



( +4 ) ( +4 



+ 4>< 



r 



+ 4 



X 



In case of the first equality, the structures in the last line on the r.h.s. can either lead 
to a permutation of field flavours or to flavour contractions. This depends on how the 
scalar lines are connected to the composite operator in the full diagram. In the other 
equalities, only the respective last term on the r.h.s. is a non-trivial permutation. 



C Permutation structures 

Since the operators in the SU (2) x SU (2) subsector are free of subtraces, the correspond- 
ing dilatation operator is given by an expansion in terms of the permutation structures 
which are defined in (12.31) . 

The definition is similar as in the M = ^ SYM case, [3], but each permutation in the 
product permutes neighbouring fields at either odd or even sites. Also, the permutation 
structures are shifted by two sides when inserted along the chain and not by one site as 
in the A/" = 4 SYM case. Permutations at even and odd sites commute with each other, 
and only do not commute if they have one position in common, i.e. if the integers in the 
argument lists of the permutation structures differ by two. The rules for manipulation 
which have to be modified compared to the ones in the A/" = 4 SYM case are 

{...,a,6, ...} = {...,6,a, ...} , |a - 6| ^ 2 , 
{a, . . . , 6} = {a + 2n, . . . , 6 + 2n} . 

We can shift all odd integers to the left and all even integers in the argument lists to 
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the right. The basis of permutation structures at four loops hence reads 



{}, {1}, {2}, {1,2}, {2,3}, {1,3}, {3,1}, {2,4}, {4,2}, (C.2) 

where the first three elements also appear at two loops. The dilatation operator decom- 
poses into three parts, acting on even, on odd and on mixed sides, respectively. These 
parts are defined by containing permutation structures with respectively only odd, only 
even and both odd and even arguments. 

D Transformations of the Feynman diagrams 

Let c(A, A){ai, . . . , a^} be the result associated with a certain Feynman graph. We can 
easily construct the contribution from the analogous graph shifted by one step along the 
chain of elementary fields from the corresponding operator. This exchanges all fields as 

A^A, Y ^Y^ , V^^V^^ coc, c*Oc*, (D.l) 

and it exchanges the colour loops of both gauge groups, i.e. it exchanges A ^ A. We 
also have to consider several sign changes due to the exchange of certain vertices and 
of the propagators of the two gauge fields. By and Vx we denote the multiplicities 
with which the corresponding propagator or vertex of type x appears in the graph. In 
X we thereby do not distinguish between the pairs of fields related as in (ID.ip . since 
the multiplicities have to consider all propagators or vertices that are related by one 
or several of the transformations in fID.ip and by reorderings of the fields. A shift by 



one side, denoted by the operator S, then transforms the contribution from a Feynman 
graph as 

S(c(A, A){ai, . . . , a^}) = [-ifA^+^A^+v^^^+yy^^^ ^(x^ A){ai + 1, . . . , + 1} . (D.2) 

Furthermore, the reflection of a given non reflection-symmetric Feynman graph con- 
tributes to the perturbation series. A reflection exchanges A A and it changes also 
the sign of all loop momenta. The type of the other field lines is kept fixed. The order 
of Y and at the composite operator and at the external lines does not change under 
reflection, if the number of interacting neighboured fields is odd. If this number is even, 
a reflection of the diagram exchanges the fields with their conjugates. If the composite 
operator is kept fixed in the latter case, a reflection of the subdiagram not involving the 
composite operator then requires a shift by one site, before it can be reattached to the 
composite operator. The so defined reflection, denoted by the operator R, transforms a 
Feynman graph as 

R(c(A, A){ai, . . . , am}) = (-l)^^^+^Ar2+v,^2+v,.,.^(^^ \){-a^, -am} , (D.3) 

and it preserves the type of the flavour permutation, i.e. its action on even or odd sites. 

The shift flD.2p by one site realizes a parity transformation. An individual Feynman 
diagram has definite parity whenever c(A, A) = ±c(A, A) with an eigenvalue given by 
this sign and the additional sign determined from its propagator and vertex content. 
There are several diagrams which individually do not have definite parity. In some cases 
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where the diagram is not reflection symmetric, the hnear combination together with the 
reflected diagram has definite parity. In cases of diagrams which involve more than three 
legs in the interaction, the diagram contains permutation structures which contribute 
to the dilatation operator at odd and at even sites (and also to the mixed component). 
Then, also the shifted diagram has to be considered, and their sum has again definite 
parity. 



E Two-loop self-energy of the scalar field 

The two-loop self-energy of the scalar fields is given by the following sum of diagrams 



+ -^ + ^^ + -^ + -!^ + '^ + -^ + ^^ + ^^ 

(E.l) 

Using the effective Feynman rules (]A.9I) and ( lA.lOl) . considering that another internal 
flavour trace arises in the first diagram giving an additional factor of 4, and that the 
additional scalar or fermion propagator respectively yield an extra factor of —ioii, the 
individual contributions read 



nf^I^ = 12i<^MiV-^- = !^(? + 4(3 - 7 + lii4ir) 



(E.2) 



The relevant two-loop integral is given in terms of G- functions (]H.3P as (72(1, 1)G(— A, 1). 
The contributions which contain quartic as well as cubic vertices all vanish according to 
the effective Feynman rules with gauge fields (lA.Qp . The same holds for the diagrams 
in the second line which contain a fermion bubble. The contributions with only cubic 
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vertices become 



AA / 2 IT- 4 , 



1 

3e 



7T 

T 



.XX 



(^3-7 + ln47r^ 



15(3 

3 ^ 



7 + ln47r) 



(E.3) 

where the integrals have first been expressed as products of two G-functions, or have 
been taken from the tables of integrals in appendix lJ.il The expression is then expanded 
in powers of e. We thereby have to keep the pole and finite parts. 

Summing up the above diagrams and their reflected copies as they appear in (IE. II) . 
the amputated scalar self-energy contribution apart from the Wick rotation (the result 
has to be multiplied by a factor P = —1) then becomes 



-•Y 



XX/ 3 3 
Tv2e ~ 2 



TT 
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y 

.2 



7 + ln47r 



— - — + -(3 - 7 + ln47r' 
4e 4 2^ ' 



(E.4) 



-l + 2e 



where the last propagator factor on the r.h.s. captures the momentum dependence. Its 
weight label indicates the exponent of ^, where p is the external momentum. The pole 
part of the above result coincides with the result in [T7] . 



F Two-loop renormalization of the six-scalar vertex 

The two-loop renormalization of the six-scalar vertex allows us to fix a sign discrepancy 
in the literatur^ parameterized by z = ±1 which affects the four- loop diagrams with 
fermion triangles, c.f. dOB]) . KTr\\ . and K28\f . 

We only compute the two-loop renormalization of the component of the scalar six- 
vertex that permutes the flavour, i.e. which has the flavour structure given by the last 



''The antisymmetric parts in the product of two 7 matrices (|A.1|) in [31] and [7] differ by a sign. 
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term on the r.h.s. of (]B.2|) . and denote it by {1}. The relevant graphs are given by 

i 
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--^MiVK(J22ij)8{l} 
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.AAl 



i (47r) 



2 fc2 



-M2K(/22a){1} 



"T4^^^> 
K(Sy){l} = ^(6AA + (A-Ar)l{l}, 



2 (47r) 



(F.l) 



— z 



4 A;2 
i (47r)2 



(M - iV)MK(/22A)(-8{l}) = -^-(A - A)A-{1} 



4 A;2 



M'K{h,,^,m-8{l})=^z^-^{l} 



where the relevant two-loop integrals are explicitly given in (lJ.2p . (IJ.3P and (IJ.lSp . The 
part which contributes to flavour permutations is determined using (lB.2p and (IB.Sp . 

Neglecting a factor = — 1 for the Wick rotation, the renormalization of the six- 
scalar vertex is then found as 



* 



^ ^ 



+ 3( ix + K 



y- 



3|l^ + >^ 



(F.2) 



= AA-z(;^-l){l}. 

Non-renormalization of the coupling and hence superconformal invariance requires z = 1, 
which corresponds to the sign choice in [3T] . 
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G Cancellation of double poles 



In the expansion of the logarithm ( 13.4p the four-loop contribution Za, is combined with 
the square of the two-loop contribution at four-loop order. This combination ensures 
the cancellation of double poles in Za, which are the remnants from two-loop subdiver- 
gences. Here, we explicitly demonstrate the cancellation of these double poles. 
The two-loop graphs are given by 




H 



MiVK(/2)( - \{} + {1}) = - ^{} + {1} 



^^M^iK(/...){} = ^l{}, (G.l) 



^K(Sk){} = -( — - + i^^l ){} 



where here the arrows indicate that we have neglected contributions to flavour traces 
but we kept the contribution to the identity. The integrals are explicitly given in (1J.2P 



f IJ.Sp . and the pole part of the scalar self-energy contribution Ey can be found in fl4.18p . 
We have also included a sign due to the double Wick rotation. The flavour permutation 
structure of the six-scalar vertex is given in (IB. 21) . In the second diagram, only the second 
substructure of fl4.10p gives a logarithmically divergent contribution to the identity in 
flavour space. 

The above results describe the sum of diagrams, in which the flrst leg is flxed to an 
odd side, c.f. the deflnition of the permutation structures in ( 12. 3p . For the flnal result, 
we have to include also the sum over even side legs, which is easily obtained employing 
the shift operator S in (lD.2p . 



The negative of the sum of the diagrams at odd and even sites with a factor \ 
for the scalar self-energy contributions then yields the two-loop contribution to the 
renormalization constant as 



i=l i i i 



where we have indicated the sum over the sites explicitly. The above result is consistent 
with [17]. It has an obvious decomposition into two parts acting exclusively on even and 
on odd sites, respectively. The square of the above result can be decomposed as follows 

-Z2 = 222,dc + ^22,5 + ^22,v + ^22,s • (G.3) 
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The individual terms are given by 



'22,dc 



j>i+3 



j>i+3 




j>i+3 




2 k 
(AA)2 1 



16 2^2 



M'^N^ K(72)'({1, 3} + {3, 1} + 2{1, 2} - 5{1}) 
({1,3} + {3,1} + 2{1,2}-5{1}) , 



i=l 



\ 




+ 



+ 



>K.>K 



^ + Ar2)M7VK(/2) K{h2A) = ^(2AA + (A - A)^)^{1} , 



22,s 



+ 




3 (4 
2~P 



-MiVi K(/2)K(Ey){l} 



A A 3 
.-(6AA + (A -Af )—{!}, 



(G.4) 

where the arrows indicate that we have neglected all contributions to the identity in 
flavour space, and also that we have only written half of all contributions, namely those 
with permutation structures with an odd entry as their first entry. It is sufficient to 
show the cancellation of the double poles for the displayed contributions only. The 
cancellations work analogously for the terms involving permutations with an even entry 
as their first entry. The gray hne in the middle of the above diagrams reminds us that 
the diagrams are direct products of the parts above and below that hne. The loop 
integrals are direct products of two loop integrals. The can be read off by contracting 
the bold gray and black line each to a point, keeping the lines which enter the bold 
gray line from below as separate external lines of the second integral. We have chosen 
the above graphical presentation to make an identification with the individual four-loop 
diagrams more evident. 

The contribution 2^22,dc cancels the double poles from those diagrams which become 
disconnected if the composite operator is removed. The corresponding four-loop contri- 
butions arc themselves squares of two-loop diagrams. They hence only contribute to the 
double poles and have not been considered in our calculation in the first place. Their 
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cancellation in combination with 222,dc is obvious. The sum of the remaining diagrams 
is given by 



^22,5 + 222.V + 222.S (^({1, 3} + {3, 1}) + {l, 2} " ^{1} + ^'^{1}) • 

(G.5) 

This result precisely cancels the double poles of the respective terms in ^4 in (15.11) in 
the expansion of the logarithm fl3.4p . 

The above cancellation works for z = 1, since according to flF.2p in this case the 
six-scalar vertex is not renormalized. In a theory where this vertex is renormalized, one 
has to remember that it is the renormalized six vertex that appears in There is 
hence an additional contribution at four loops which comes from the two-loop part of 
the respective vertex renormalization constant in Z2 itself. If this is taken into account, 
double pole cancellation should work also in this case. 



H C-functions 

To compute the required integrals in dimensional reduction in D dimensions, we make 
use of G-functions. The scalar G-function is defined as follows 

^^^^^ . , d^/ 1 r(a + /3-f)r(f -«)r(f 



i2n)D J - p)2/3 (4^)f T{a)T{P)T{D - a - /3) 

(H.l) 

where the integral is defined in Euclidean space. To analytically continue to Minkowski 
space, we have to identify real time as t = —ir in terms of Euclidean time, and hence 
multiply each loop integral by i. The non-trivial tensor G- functions can be obtained 
from this by introducing traceless symmetric products of the loop momentum in the 
numerator. Traceless symmetric tensor indices are embraced by parentheses. We then 
define 

p2(.+,-f) ^^ = J2^J P^ii-py^ ,2=1 • ^^-2) 

For a one-loop integral with a single or two contracted momenta in the numerator we 
define the G-functions without parentheses as 




jj-G(a,/3) 



-<^= o, G,{a,P) , Gi(a,/3) = ^(G(a,/3)-G(a,/3-l) + G(a-l,/3)) , 

-- — -J--^G2(«, (3) , G2(«, P) = Gi(a, /3) - G{a -l,f3), 

(H.3) 

where an arrow on a propagator denotes the appearance of the momentum running along 
this line in the numerator, and two arrows of identical type denote a contraction of the 
respective momenta. The function Gi is identical to the function G(i). The function 
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for a traceless symmetric product of two equal momentum factors in the numerator 
explicitly reads 



G(2)(a,/3) 



D 



G{a, /3) + G{a -2,(3) + G{a, (3-2) 



+ 2 



D-2 
D 



G{a -1,(3)- G{a, (3 - 1) - G{a - 1, (3 - I) 



(H.4) 



I Triangle rules 

Defining the functions 

A{a,(3) 
Cia,P) 



aG{a + l,(3) + (3G{a,(3 + l) 



a + (3 + 2-D 



a 



a+P+2-D 
the scalar triangle rule is given by [38|l39] 




A{a,(3) 



+ C{a,(3) 



a + l 




+ C{(3,a) 




13+1 



Defining the functions 



A±{a,(3) = A,{a,(3)±A{a,(3) 

^ , ^. [a- (3)G[a,f3)-aG[a + \,f3-\)+ i3G[a-\,^ + \) 
Ax[(x,(3) , 

A{a,(3) 
C{a,P) 



2{a + (3 + l-D) 

aG{a + l,(3) + PG{a,(3 + l) 

2{a + (3 + l-D) ' 
a 



a + (3 + l-D ' 

the triangle rule with a single momentum in the numerator reads [39] 




A+(a,/3) 



a + l3 + l-^ 



A_(a,/3) 



a+/3 + l-^ 



+ C{a,(3) 




C{P,a) 




/3 + 1 



(I.l) 



(1.2) 



(1.3) 



(1.4) 
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We then immediately find from (]I.2|) and f ll.4|) also the following rule 





+ 



(A,(a,/3)-iA(a,/3))"-;^^ 
A(a,/3) 



1 



2(a + (3 + l-D) 



a + + 



+ 



Q+/3 + l-4j 



+ (C'(a,/3)-C(a,/3)) 




+ (C'(/3,a)-C(/3,a)) 




/3+1 



+ C(a,/3) 



Q + l 





Defining the functions 

A2(a,/3) = 



1 



2(a + /3 + 2-L)' 



1 



■ {a{G{a + 1,13 -1) -G{a + 1, /?)) 

+ f3{G{a-l,f3 + l)-G{a,f3 + l)) 
+ {a + f3-2)G{a,f3)) , 



2(a + /3 + 2-D) ' 

we obtain the triangle rule with two contracted momenta in the numerator as 




D{a,(3) 



+ G{a,f3) 




+ 




13-1 






G{f3,a) 




/3 + 1 



(1.5) 



(1.6) 



(1.7) 



J Tables of integrals 

In the following we list all required integrals. The integrals are regularized by dimen- 
sional regularization in Euclidean space with dimension 



D = 2(A + 1) 



(j.i) 



The parameter A in this appendix always assumes the above value. It has nothing to 
do with the 't Hooft coupling that appears in the main text and which is also denoted 
by A. The integrals have a simple dependence on the external momentum p^. 
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J.l Two- loop integrals 

The simplest two-loop integral is given by 

Ha) = = ^(1, l)G{l - A, a) . 



(J.2) 



It is proportional to 2(}-2\) ■ the main text we abbreviate I2 = -^^2(1)) which is the 
simplest logarithmicafly divergent integral in three dimensions. 

The two-loop integrals with a bubble and two contracted momenta in their numer- 
ators read 

h2A= lQ\ =G2(1,1)G(1-A,1) , 




I22B 
I22C 
I22D 






Gi(l,l)G(l-A,l) , 
Gi(l,l)G2(2-A,l) 
G'(1,1)G2(2-A,1) 



(J.3) 



They are proportional to ^2{i-2X) ■ 

We also need the following 2-loop scalar integral which is proportional to 2(a+2-2A) 



/20a (a) 



G(l, l)(A(a, 1) + C{a, l)G{2 - A, a + 1) 
+ C(l,a)G(a + l-A,2)) , 



(J.4) 



where one propagator has weight a, i.e. its squared momentum in the denominator has 
exponent a. 

More difficult to obtain are the two-loop integrals in which the central line has a 
propagator weight differing from one. To compute the needed graphs, we work with p- 
space weights in p-space, not with p-space weights in x-space as in the rest of the paper. 
The corresponding special 2-loop graph in p-space looks the same as in x-space, but two 
lines which are connected at one point are fused with the G-functions, while the weights 
of two lines forming a loop can be directly added to yield the weight of the single line 
which replaces the loop. First, we can use the point transformations described in [ID] 
to relate 




2A+1-Q-/3 




2A-/9 



^( + 1-A 



(J.5) 



where in the ffist, second and third equality the conformal transformation of inversion, 
the adding of an external line to make the vertex unique, and the conformal transfor- 
mation of inversion, have been respectively used. We can then use the transformation 




a + l3-\ 
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to evaluate the following integrals 



A A 

A = Q =Wl) = ^(-^-4(l-, + l„4.)). 

The first integral is finite, and thus its finite part is identical to the finite part of one of 
the required integrals. We obtain 



1 .2 



20a = ^ = (J.7) 



up to terms of order 0{e). 

The second of the above integrals is correctly reproduced by the F-function of [H] as 
3A — 2, 1), which is given as a converging double series for the given values. We 
can then use the F-function to also compute another required integral with a different 
^-dependence of the central hne. It reads 



(47r; 



i,op= (-^ =y^F(l,-A,l) = -i-^(^ + 4(l-7 + ln47r)) . (J.8) 



The above results are used to compute the following integrals which are proportional 

to 1 

'j'-' p2(2-3A) 

1/1 7r2 , X 

+ 6- — -27 + 21n47rj , 



(87r)ne 2 
h2,= M =^/2o/3 = T^(^ + 2-27 + 21n47r 



(J.9) 



/22,= M =^(2G(l-A,l)G(2-2A,l)-/2o^) = -l-^4 



1 /. 

ttPV" 2 



'225 — t 1 — -^22a — -^22/3 — t:"^ I 4 — 
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W = = 2(^22/3 - ^2(2 - A, l)G{2 - 2A, 1) - 6-1(2 - A, 1)G'(2 - 2A, 1)) , 



T — 13^ — T 



222/3€ ~ -^22/3 • 



(J.IO) 

We need the following integrals with a contraction of one loop momentum in the 
numerator with the external momentum which are proportional to p2(«+i-2A) 



hla{0L) 



l)G{a + 1 - A, 1) + hoaia) - hoaia - 1)) , 



l)G(a + 1 - A, 1) - hoaia) - hoaia - 1)) , 



/2ic(«)= =Gi(l,l)(A_(«,l) + A+(aa) + C(a,l)G'i(2-A,a + l) 

-C'(l,a)Gi(a + l-A,2)) , 

/2w(a)= (_) =^(-G'(l,l)G'(2-A,a) + G'(l,l)G'(a + l-A,l) + /20a(«)) , 



/2ie(«)= M =-(-G(l,l)G'(2-A,a) + G'(l,l)G'(a + l-A,l)-/2o„(a)) , 



(J.ll) 

We need the following integrals with two contracted loop momenta in the numerator 
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which are proportional to 



(a+l-2A) 



l22a{a)= (-^ =i(G'(Q;,l)G'(Q; + l-A,l)- 720a(Q;) + /20a(Q;-l)) , 



I22M= (-^ =i(G'(l,l)G'(a,l) + G'(a,l)G'(a + l-A,l) 

-G'(1,1)G(q; + 1- A,l)) , 

722c(a) = (-^" = ^(-^(1, 1) + l)G'(a + 1 - A, 1) 

+ (?(!, l)G'(a + l-A,l)) , 

^22d(a) = — \ = h2a{a) - l22b{a) , 



722e(a) = = 1{G{1, l)G{a, 1) - G(l, 1)G(2 - A, a) + /20a(a - 1)) , 



(J.12) 



-^22/ (tt) = L ] = /22a (") - -^22e(") , 



l22g{a)= {—] =^{-G{l,l)G{a,l) + G{l,l)G{2-X,a) + l20a{a-l)) , 



l22h{a)= = l{-G{l,l)G{a,l) + G{a,l)G{a + l-X,l) 



2 



G{l,l)G{a + l-X,l)) , 



l22i{a) = = li-G{l, l)G{a, 1) - G{1, 1)G{2 -X,a) + hoaia - 1)) , 



/22,(a)= (— ) =^(G(l,l)G(2-A,a) + G(l,l)G'(a + l-A,l)-/20a(«)) , 



We need the following integrals in which two loop momenta in the numerator are 
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contracted with the external momentum and which are proportional to 2{a-2x) 



-^2110(0;) = 




-(Gi(l, a)Gi{a + 1 - A, 1) + G{1, a)Gi{l, a + l-X) 

- hlaipi) - hlaia - 1)) , 



-^2116(0;) 



;(G'i(l, a)Gi{a + 1 - A, 1) + hida) - hidot - 1)) , 



hiic{oi) = 



l)Gi(2 - A, a) + ^(1, 1)6-1(1, a + 1 - A) + 72ia(«)) , 



;-G(l, l)Gi(2 - A, a) + G{1, 1)^1(1, a + 1 - A) - hia{c^)) , 



-(Gi(l, a)Gi{a + 1 - A, 1) - hUa) - hiM - 1)) , 



/2ii/(a) 



If 1 



(G(l, l)G'i(a, 2 - A) - G(l, l)G^{a + 1 - A, 1) + hi,{a)) , 



-(G(l, l)G^{a, 2 - A) - G{1, l)G^{a + 1 - A, 1) - /2ib(a)) , 



-{-G,{1, 1)G,{2 - A, a) - Gi(l, l)Gi{a + 1 - A, 1) + hida)) , 



[-G,{1, 1)G,{2 - A, a) - G,{1, l)G,{a + 1 - A, 1) - hu{a)) , 



-{-G^{l, l)Gi(2 -X,a) + G^{l, l)G,{a + 1 - A, 1) 
+ G{1, l)Gi(l, a + 1 - A) - hidia)) , 



(J.13) 

We need the integrals in which two loop momenta in the numerator are contracted 
with each other and one loop momentum is contracted with the external momentum. 



41 



The integrals are proportional to 2(a-2x) and read 



h21ac{0l) 



-(^(l, a)Gi{a + 1 - A, 1) - hida) + hida - 1)) , 



l221ad{0l) 



-{G{a, l)Gi(l, a + 1 - A) - hui^) + hui^ - 1)) , 



l221ae{oc) 



-{-G{a, l)Gi{a + 1 - A, 1) - huia) + hiei^ - 1)) , 



-^22166 (a) 




-(-^(1, l)G'i(a, 1) - ^(a, 1)^1(0; + 1 - A, 1) 
+ G'(l,l)G'i(a + l-A,l)) , 



-^2216c(a) 



-(G(l,l)Gi(l,«)-Gi(l,l)G(a,l) 

+ ^1(1, a)G'i(a + 1 - A, 1) + ^(l, l)Gi{a + 1 - A, 1)) , 



-^2216d(a) 



-^22166 («) 



^(^(l, l)G{a, 1) + G{a, 1)^1(1, « + 1 - A) 
- ^(1, l)G'i(l, «+ 1 - A) - 6-1(1, l)G'i(a + 1 - A, 1)) , 

^(-^(l, l)G'(a, 1) - l)G'i(a + 1 - A, 1) 
+ G'i(l,l)G'i(a; + l-A,l)) , 



l22lcc{0i) 



-(-^(1, 1)^1(1, a) + 6-1(1, l)G{a, 1) + 61(1, a)Gda + 1 - A, 1) 
-G'i(l,l)G'i(a + l-A,l)) , 



h2\ce{0L) 



-{Gx{l, \)G{a, 1) - G{(x, \)GdoL + 1 - A, 1) 
- 6-1(1, 1)6-1(0; + 1- A, 1)) , 



h2\dc{0l) 



h2\ac{oi) - l22\hc{oi) , 



h2\M{pt) 



h21ad{0L) - l221bd{a) , 



(J.14) 
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l221ebM 



-G{1, l)G,{a, 1) + G{1, 2 - A) - hib{a - 1)) , 



-(Gi(l, l)G{a, 1) - 6-1(1, 1)61(2 - A, a) + h^a - 1)) , 



h21ee{0i) 



'-Gi{l, l)G{a, 1) + Gi(l, l)Gi(2 - A, a) 
+ G(l,l)Gi(a,2-A) + /2ie(a-l)) , 



l221fe{a) 



h21ae{a) - /221ee(a) , 



h21gd{a) 



(-Gi(l, l)G{a, 1) + Gi(l, l)Gi(2 -X,a) + huia - 1)) , 



'221fea 



-(-G(l, l)Gi(l, a) + Gi(l, a)Gi{a + 1 - A, 1) 

+ G{a, 1)61(1, q; + 1 - A) - 6(1, 1)61(1, a + 1 - A)) , 



where e.g. we have the relations 



h21ad(o:) - l22lae(o:) = /22a(a) , 
l221bd{a) - /2216e(a) = -^226(q:) , 
l221ed{a) - /221ee(a) = -^22e(a) , 



(J.15) 



(J.16) 



We need the integrals in which four loop momenta in the numerator are pairwise 
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contracted with each other. The integrals are proportional to 2(a-2x) and read 



-^222oe(a) 



-(-(7(1, l)G2{a, 1) + l22aia - 1) + ^(1, 1)^2(2 - A, a)) , 



h22af{0l) 



-G,{a, l)G'2(a + 1 - A, 1) - 722/(a) + /22/(« - 1)) , 



-^222aft(Q;) 



-(^(1, l)G2{a, 1) - G2(a, l)G(a - A, 1) 

- Gi(a, l)G2{a + 1 - A, 1) + ^(1, l)G2{a + 1 - A, 1)) , 



-^222ai (a) 



-(^(l, l)G2{a, 1) + /22a(« - 1) + ^(1, 1)G2(2 - A, a)) 



' 222ai 



[a) 



-{-G{a, l)G2{a + 1 - A, 1) - l22j{a) + /22,(« - 1)) , 



-^222b/(tt) 



/22265(a) 





1 



-Gi(l, l)G'i(a, 1) - Gi{a, l)G2{a + 1 - A, 1) 

+ G(l, l)G2(a + 1 - A, 1) + Gi(l, l)G{a - A, 1)) , 



= -(^1(1, l)G'i(«, 1) + Gi{a, l)G{a - A, 1) 



-Gi(l,l)G(a-A,l)) , 



-^222&/i(a) 



^(^(l, l)G'i(l, a) + G,{1, a)G2ia + 1 - A, 1) 
+ Gi(1,1)G'2(q; + 1- A,l) , 



l222bi{0i) 



= 1^{G^{1, l)G'i(l, a) + 7226(a; - 1) - ^(l, 1)G{1 - A, a)) , 



l222bj{a) 



-(-6-2(1, 1)6(0;, 1) - G{a, 1)62(0; + 1 - A, 1) 

+ 62(1, 1)6(0; - A, 1) + 61(1, 1)62(0; + 1 - A, 1)) , 



(J.17) 
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-(-6-1(1, 1)6-1(1, a) + Gi(l, a)G2{a + 1 - A, 1) 
+ G-i(l,l)G-2(aH-l-A,l)) , 

^(Gi(l, l)G-i(l, a) + /22c(a - 1) - Gi(l, 1)G-(1 - A, a)) , 

^(-G-i(l, l)G-i(l, a) + l22d{a - 1) + G-i(l, 1)G-(1 - A, a) 
+ ^(1, 1)6-2(2 -A, a)) , 

l222af — -^2226/ , 

^(G-i(l, l)G-i(l, a) - G-i(l, a)G{a - A, 1) 

- G{a, l)G2{a + 1 - A, 1) + 1)6-2(0; + 1 - A, 1)) , 

^(Gi(l, l)G-i(l, a) + h2d{a - 1) + (7(1, 1)^2(2 - A, a) 
+ G-i(1,1)G-(1-A,q;)) , 

^(-G-i(l, l)G-i(o;, 1) + G-i(l, l)G-2(2 - A, a) + 722/(q; - 1)) , 

^(G(l, l)G-2(a, 1) - G-i(l, 1)G-i(q;, 1) + 722e(a - 1) 
-Gi(l,l)G-2(2-A,a)) , 

i(-G2(l, l)G-(a, 1) + /22i(a - 1) + ^2(1, 1)G(1 - A, a) 
+ G-i(l,l)G2(2-A,a)) , 

^(Gi(l, l)Gi(a, 1) + /22/(a - 1) + Gi(l, 1)6-2(2 - A, a)) , 

(J.18) 
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J. 2 Three-loop integrals 

J. 2.1 Three- loop integrals with central cubic vertex 

The positions for the numerator momenta are indicated as foUows 




(J.19) 



The direction of the numerator momenta (arrows) is such that the label is on the r.h.s. 
if one follows the momentum. The integral with two loop momenta in the numerator 
which are contracted with the external momentum reads 




= A_(l, l)(/2ii;,(2 - A) - /2ii/,(2 - A)) + A+(l, l)/2m(2 - A) 

+ C{1, l)(/3iic(2) + G{2, l)/2iu(2 - A) - 6-1(1, 2)/2iie(2 - A)) 



1 



(Stt) 



(J.20) 

We need one integral in which two loop momenta in the numerator are contracted with 
each other and one loop momentum is contracted with the external momentum. It reads 



'321cci/e 




[ - G{1, l)72ic(l - A) - (7(1, 1)^1(2 - A, 1)^(3 - 2A, 1)] 



2 

(Stt) 



1 / Tl\ 



(J.21) 
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The integrals in which four loop momenta in the numerator are pairwise contracted read 




h22oacde= UU\ = - [ - Gi (1, 1)2G2(2 - A, 2 - A) - 2^ (1, l)/22a(l " A)] 



h22cadce — ( > ^ < 




h22cadde — [» ' < 




1 TT 



(87r)= 



^( -- - l + 37-ln2567r^) , 

2 \ £ / 



- G{1, 1)7226(1 - A) + G{1, 1)G'2(2 - A, 1)G{2 - 2A, 1) 

+ ^G(1,1)2G(2-2A,1) 



l^.(-i-4 + ^ + 37-ln 256.3), 



- [ - G{1, 1)G'2(1, 1)^(2 - A, 1 - A) + 1)722,(1 - A) 
-G(l, 1)7226(1 -A)] 

1 ""(-l-S+y + 37-ln2567r3) , 



(87r)3 2 



(J.22) 
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'322cac(ig 




= -(Ai(l,l)-iA(l,l))7222ai(2-A) 



^ -A(l, l)(/222<ii(2 - A) - h22M{2 - A)) 



2(3 -L>) 

+ ((7(1, 1) - C(l, 1)) ( - ^(1 - 2A, 2)l222dhil) + Gii2 - 2A, 2)722idc(l) 

+ G(2)(l, 2) (/222e/(2 - A) - -^/22^(1 - A)) 
-G'i(l,2) 7222e,(2-A)) 

+ C(l, 1) (^2(2 - 2A, 2)/22im(1) - 6-1(2 - 2A, 2)(/22i66(l) - /221e;,(l)) 

- G'(2)(3 - 2A, 2) (72ih(1) - -^/22i(l)) - ^G{2 - 2A, 2)722^(1) 

+ G'(2)(l, 2) (7222a/(2 - A) - -^722c(l - A)) - ^(l, 2) 7222a,(2 - A)) 



1 / 2 9 

rTT -6 + -TT^ 



'322cascd — t"^ 




= -(Ai(l,l)-iA(l,l))7222<ie(2-A) 



^ -A(l, 1)(7222*(2 - A) + Gi(l, l)Gi(l, 2 - A)) 



2(3-L>)' 

+ ((7(1, 1) - C(l, 1)) (g{1 - 2A, 2 ) 72226^(1) + G,(2 - 2A, 2)722i„e(l) 

+ G'(2)(l, 2) (7222e/(2 - A) - -^722/^(1 - A)) 

-(71(1,2 ) 7222^(2 -A)) 

+ C(l, 1) ( - G2(2 - 2A, 2)72216d(l) - ^1(2 - 2A, 2)(722m(l) - 7221e6(l)) 

- G^2){S - 2A, 2) (7211,(1) - -^722,(1)) - ^G{2 - 2A, 2)722,(1) 

+ G(2)(l, 2) (7222a/(2 - A) - -^722c(l - A)) - ^(l, 2) 7222<i/(2 - A)) 

(J.23) 
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'322c6dce 




Iz22cbecd — 




'322cc/de 




'322c(iede 




'322cdee/ — 




'322c(i/eg 




^ [6-1(1, 1)/22,(1 - A) + 6(1, 1)/226(1 - A) 

-G(1,1)G2(1,1)G(2-A,1-A)] 



- [ - 61(1, l)/22a(l - A) + 6(1, 1)7226(1 - A) 

- 61(1, 1)2G'2(2 - A, 2 - A) - 61(1, l)722e(l - A)] 
(Stt^ 



- 1 _ ^\ 



■ [ - G{1, 1)^1(1, 1)6-2(2 - A, 2 - A) - Gi(l, l)/22/(l - A) 

-G(l,l)/22d(l-A)] 



(Stt)^ 



1 



[61(1, 1)^62(2 - A, 2 - A) + 26-1(1, l)722e(l - A) 
- 26(1, 1)/22.(1 - A)] 
L_|(l + 9_3, + ln256.3), 



(Stt) 



- [ - 6(1, 1)61(1, 1)62(2 - A, 2 - A) + 6(1, l)/22i(l - A) 

+ 6(1, l)l22h{l - A) - 61(1, 1)722/(1 - A)] 

1 / vr^ 

TT — 1 

(87r)3l 4 



^6(1, 1) [722.(1 - A) - ^6(1, 1)6(2 - 2A, 1) 



+ 61(2- A, 1)6(2 -2A,1) 



(8' , 

(J.24) 

We can then evaluate certain integrals as linear combinations of the above integrals. 
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They read 




(J.25) 

J. 2. 2 Three- loop integrals with central triangle 

The positions for the numerator momenta are indicated as follows 




(J.26) 



The direction of the numerator momenta (arrows) is such that the label is on the r.h.s. 
if one follows the momentum. 

The integrals in which two loop momenta in the numerator are contracted with each 
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other read 



h2tad — ^' 



1) [ - G(l, 1)G(2 - A, 1) + /20a(l - A) + hoa] = T^y ' 



'32ta/ 



'32tce 



= A(l, l)/22^i + 2C(1, 1)(G'(2, l)/2id(2 - A) + ^1(1, 2)/22/(2 - A)) 



-h2taf — l)/20a + 2/32toci 





(87r)= 



-47r) 



(J.27) 

We also need one integral, in which two loop momenta in the numerator are con- 
tracted with each other and one loop momentum is contracted with the external mo- 
mentum. It reads 



'321tces 



- [ - 1)G2(2 - A, 1)G(A, 2 - 2A, 1) - Gi(A, l)722e(l - A) 

~ -^32tce] 

1 

r27r . 



(Stt) 



(J.28) 

We also need several integrals in which four loop momenta in the numerator are 
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pairwise contracted. They read 



' 322tadce 




'322taddff 




■ l)2Gi(2 - A, 1) - ^(1, 1)G2(1, 1)G{1, 1 - A) 

- 6-1(1, l)/22e(l - A) + 6-2(1, l)/20^] 



(87r)= 



TT 



- [G{1, 1)61(1, 1)62(2 - A, 1) + 6(1, l)722d(l - A) + 61(1, 1)722/3 
+ 6(1, l)/225] 

0{e), 



'322taecd 



I 322ta fdg 




- [61(1, 1)^61(2 - A, 1) - 61(1, l)/22e(l - A) 

+ 6(l,l)/22fe(l-A)-6i(l,l)/22/3] 
1 /5 TT^ 



-A_(l, 1) 722275 + A+(l, l)7222aT 

- 6(1, 1) ( - 6(2, l)722idd(2 - A) - 61(1, 2) 7222d/(2 - A) 
-6i(l,2)722w(2-A) 

- 6(2)(1, 2) (7222a/(2 - A) - -^722c(l - A))) 



(Stt)^ 



r57r 1 



12/ ' 



(J.29) 
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A2(l, l)/225 

-L'(l,l)(-G'(l,l)G'i(l,l)G'i(l,2-A) 
-G(l,l)Gi(l,l)Gi(2-A,l) 
+ Gi(l, 1)/225(1 - A) - ^(1, l)/22<i(l - A) 

-G'i(l,l)/22/(l-A)) 
+ C(l, 1) (G'(2)(2, 1) (7222e,(2 - A) - 1/22.(1 - A)) 

+ -^^(1, 1)/22,(1 - A) - 6-1(2, 1)W(2 - A) 

+ G'(2)(2, l)(/222e/(2 - A) - -^122/^(1 " A)) 

+ ^G{1, 1)/22^(1 - A) - Gi(2, l)/222e/(2 - A))) 



i [G(l, l)Gi(l, l)Gi(2 - A, 1) + G(l, l)/225 



(Stt) 



+ 61(1,1)722/(1 -A)] 



A2(l,l)/22-y 

- 2D(1, 1)(G'(1, 1)6-1(1, l)Gi(l, 2 - A) - Gi(l, l)/22g(l - A)) 

+ 2C(1, 1) (G(2)(2, 1) (7222e,(2 - A) - 1/22.(1 - A)) 

+ 1g(1, 1)/22.(1 - A) - 61(2, l)/222e,(2 - A)) 



^ :-27r), 



(87r)= 

(J. 30) 
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^ [ - 1)^2(2 - A, 1) + (7(1, 1)^(1, l)Gi(2 - A, 1) 

-Ca, 1)722,(1 - A) +^1(1, 1)722/3] 

_i_^(i + 2-3, + ln256.3), 



^ [^(1, 1)G2(1, 1)^(1 - A, 1) + ^(1, l)722b(l - A) + 6-1(1, 1)722,3] 



^7rf- + 5-^-37 + ln2567r=^') , 



^ [G(l, 1)61(1, 1)62(2 - A, 1) - 6(1, 1)722„ + 61(1, 1)722/3 



+ 6(1,1)722^(1 -A) 
TT r 2 \ £ / 



(Stt) 



i [ - 6(1, 1)61(1, 1)62(2 - A, 1) + 6(1, 1)722,(1 - A) + 61(1, 1)722/3] 
1 n (\ 



-(- - 2-37 + ln2567r=^) , 

2 V £^ / 



\ [6(1, 1)^62(2 - A, 1) - 6(1, 1)722. + 6(1, 1)722,(1 - A)] 
^.(-l-4 + !^ + 3,-ln256.3), 



i [6(1, 1)62(1, 1)6(1 - A, 1) - 61(1, 1)^61(2 - A, 1) 



+ ^2(1. 1)720/3 + 61(1, l)722e(l- A)] 



^ [ - 61(1, 1)^61(2 - A, 1) + 61(1, l)722e(l - A) - 61(1, 1)722^] 

(8^(~ 2) ■ 

(J.31) 
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^ [6-1(1, 1)^61(2 - A, 1) - 6-1(1, l)/22/3 - G(l, l)/22e(l - A)] 



1 TT 
^2 



- ^ - 1 + 37 - In 2567r^) , 



^ [ - 6(1, 1)6-2(1, 1)6(2 - A, 1) + 6(1, l)/22a + 62(1, l)/20a(-A) 

-6i(l,l)/22a(l-A)] 
1 TT/l 



B7r)3 2 



+ 4 - y - 37 + ln2567r=^) , 



^ [61(1, 1)^61(2 - A, 1) + 61(1, 1)722/3 + 6(1, 1)7225 

+ 6i(l,l)722a(l-A)] 
1 TT/l 



B7r)3 2 



+ 7 - y - 37 + In2567r3) , 



i [6(1, 1)62(1, 1)6(2 - A, 1) + 6(1, 1)722. - 6(1, l)72o„(-A) 

+ 6i(l,l)722a(l-A)] 

(8^f(i + ^-T-3^ + ^^256.3), 



^ [ + 61(1, 1)^61(2 - A, 1) - 61(1, 1)722^ + 6(1, l)722a 



-6i(l,l)722a(l-A)] 
1 /5 _ TT^ 



(J.32) 
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-Gi(l, l)/22/3 — ^(1, l)/225 + h22tdfdg + h22taddg ~ h22tafdg 
1 



(Stt) 



13 27r2 



' 322t cecg 



322tcedg + -^322tagce — 



1 /13 27r2 



(Stt) 



'322taecff 



h22tcecg + h22tcgde 



(Stt, 



1 / TT^ 

27r(3 



'322tc/cff 



h22tcdcg ~ 



1 TT / 1 47r^ \ 

3m».» = 7^2(j-n-3T + ^ + ln256.') , 



'322tc/dfl 



1 TT / 1 TT^ \ 

322tdfdg + h22tafdg = (-g^j3 2"(, + — + 37-ln 2567r^j , 



(J.33) 



J. 2. 3 Three- loop integrals with central square 

The positions for the numerator momenta are indicated as follows 



(J.34) 



The direction of the numerator momenta (arrows) is such that the label is on the r.h.s. 
if one follows the momentum. 

The required integral can be easily computed by promoting it first to a logarith- 
mically divergent four-loop integral by adding a scalar line, and then cancelling one of 
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its scalar lines (after shifting external momenta thereby taking care not to produce IR 
divergences). We hence find 



'322sc/de 




l322cdgef 



(8x) 



(J.35) 

where the equality between the two integrals is only valid for the finite parts, they may 
differ at 0{e). Two arrows of the same type at a single hne indicate a factor of squared 
momentum along that line in the numerator, and hence cancellation of the propagator 
(which in the above graph means shrinking it to a point). 

The required integrals with six pairwise contracted loop momenta in the numerator 
are then found as 



'3222scede4f 



'3222sc/dede 



- [^(1, 1)26-2(2 - A, 2 - A) + G{1, 1)'G'2(1, 1) - G{1, l)h2c] 

+ ^[- G(l, lfG2{2 - A, 1) + G(l, l)Gi(l, 1)G2(2 - A, 1) 
+ G{1, 1)62(1, 1)6(1 - A, 1) + Gi(l, 1)722/3] 

+ G'i(l,l)722a(l-A) 



1 /I , TT^ 
TT - + 1 + — 




37 + ln2567r3) , 




+ 



h22tcede + Iz22tcedg + -;^h22scfde 



(Stt) 



f - - 6 + TT^ - 37 + In 2567r2') 



(J.36) 
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J. 2. 4 Three-loop integrals for the nearest-neighbour interactions between 
two internal legs 

The above results for three-loop integrals with non-trivial numerators are required to 
determine the following three-loop integrals which enter fl4.15p 



'3fb 



-2G(1, 1)2^2(1,1) , 



'3gb 



6-2(1, l)/2o/3 = + 2 - 37 + In 25671^) 



'3gt 



1 / 

-6*2(1, ^)hoi3 + h22tcdde = /g^Ng Trl^S — 



'3gs 



4 [62(1, l)/20/3 ~ 2/322tccZcZe + 2I^222scededf — I?,222scfdede\ — 0{e) 



'3gn 



2 [62(1, l)l22h + 62(1, l)/22fe — Im2tcgde — -^322tcdeg + I-i22tcedg\ 
1 



vr 



'3gc 



G'l(l,l)(/222/3e ~ -^22q ~ I225 + -^222/3i) 
+ Iz22caecd — I'i22cacde — I'i22ccfde + -^322ccde/ 
1 / TT^X 



'3gv 



4 [ - 6(1, l)/222;3. + ^(1, l)/22j = 7^7r(24 - 271 



Thereby, the effective Feynman rules ( lA.QI) and OA. 101) have been used. 



(J.37) 
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J. 3 Four- loop integrals 

In this appendix we list the appearing logarithmically divergent four-loop integrals. In 
contrast to the three-loop integrals, we only need their pole parts and directly give the 
expressions in which the appearing subdivergences have been subtracted. 
The simplest scalar four-loop integrals are given by 



J. 3.1 Four-loop integrals with three bubbles and two contracted momenta 
in the numerator 

We also need the following integrals 




= K(K(G'(1, 1)^G'(1 - A, 1)G'(2 



3A,l))-K(72)/2) 




(J.38) 




K(G'(1, 1)^2(1, 1)^(1 - A, 1)^(2 - 2A, 1 - A)) 
- K(G(1, 1)G(1 - A, 1))G2(1, 1)G(1 - A, 1)) 




(J.39) 




K(G'(1,1)2G'2(1,1)G(2-3A,1)) 



1 TT 
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J. 3. 2 Four-loop integrals with two bubbles 

We need the following integrals with two bubbles and two contracted momenta in the 
numerators 




^42bbide = = K{G{1, lyGiil, 2 - 2A)G'2(3 - 3A, 1)) 

1 / TT^N 



- 



(87r)n 4s J ' 

/42bb3ad = = - l)Gi(l, 1)G2(2 - A, 1)G(2 - 2A, 1 - A) 

-K(G'i(l,l)G'2(2-A,l)/2) 
1 / l^ 




/42bb3as = ^ 1 = - ^0(1, 1)6-1(1, 1)6-2(2 - A, 2 - 2A)G(2 - 3A, 1)) 





Wscd = = - K(G(1, 1)G2(2 - A, 1)6(2 - 2A, 1)6(2 - 3A, 1) 

-K(6(l, 1)62(2- A, l))/2) 

" (8^(~ 2^2) ' 

W4a6 = = K(6(l, 1)62(1, 1)6(1 - A, 1)6(2 - 2A, 1 - A)) 

-K(62(l,l)6(l-A,l))K(/2) 
1/11- 




7r)4 V 4^2 £ 



W4ad = = - K(6i(l, 1)62(2 - A, 1)6(2 - 2A, 1)6(2 - 3A, 1) 

-K(6i(l, 1)62(2- A, l))/2) 
1 / 





(J.40) 



(87r)n 4:eV ' 

h2bb4de = = K(6(l, 1)62(2 - A, 1)6(2 - 2A, 1)6(2 - 3A, 1) 

-K(6(l, 1)62(2- A, l))/2) 
_ 1 1 
" (8^2i2 ' 

We also need the following integrals with two bubbles and four pairwise contracted 
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momenta in their numerators 




Wbsacbe = = - ^(^(1, 1)G,(1, 1)G'2(2 - A, 1)G(2 - 3A, 1) 

-K(G'i(l,l)G'2(2-A,l))72) 
1 ^-±\ 



(87r)n 4£2; ' 

Wbsaebc = = - K(G'(1, 1)^1(1, 1)^2(2 - A, 1)G{2 - 3A, 1) 

-K(G'i(l, 1)^2(2- A, 1))72) 
1 / ±^ 




(87r)n 4£2y ' 

h22bb3becd = \^ = - K(G(1, 1)G2(1, 1)^2(2 - A, 1)G{2 - 3A, 1) 

-K(G'2(1, 1)6-2(2- A, l))/2) 

1 ^-J-^ 




(87r)n 4£2y ■ 

J. 3. 3 Four-loop integrals with one bubble 

We need the following integrals with one bubble and two contracted momenta in the 
numerators 



U2had 




K{G{1, l)/22e(l - A)G(2 - 3A, 1) = ^ ( - ^) , 




W=\X7 =-K(G(l,l)/22i(l-A)G'(2-3A,l))= ^ ^ 



4 ^ ' 



h2bbe =\X\ / = K{Gi{l, l)/22e(l - A)G(2 - 3A, 1) + K(G'i(l, 1)^2(2 - A, l))/2) 

(87r)n 4£2; 

(J.42) 

The required integrals with one bubble and four pairwise contracted momenta in 
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their numerators read 



'422cb6ad&e 




'422cb6&c<ie 




'422cb66ecd 




U22ch7adbd 




K(G(1,1) 

l222ah 

{l)G{2 - 3A, 1) - K{h22ah{l))l2) 
1/1 1 . 1 ^2, 



1 TT 

2 + ¥ 



K(G(l,l)/222j,,(l)G(2-3A,l) 



^7r)n 16£/ ' 



- K(G(1, l)/222v(l)G(2 - 3A, 1) = ^ J ( - i + ^) , 



- K(G(1, 1)G(2 - 2A, 1 - A)(-G(l, 1)^2(2 - A, 1) 

+ G'i(l,l)G'i(2-A,l) 

+ G'(1,1)G2(1,1)) 
- K(G'2(1, 1)^(1, 1) + (72(1, 1)^(1 - A, 1) 

+ G'i(l,l)G'2(2-A,l))72) 



1 1/1 
4^^ £V2 ^"S" 



(J.43) 



J. 3. 4 Four-loop integrals with central quartic vertex 

The positions for the numerator momenta are indicated as follows 





a \. 


I d 






b ) 


\ c 





(J.44) 



The direction of the numerator momenta (arrows) is such that the label is on the r.h.s. 
if one follows the momentum. 

With the external momentum entering the central vertex and exiting the lower ver- 
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tex, the required integrals become 



'422qASad 




G(2 — 3A, l)l322tadce 



11/ 1 



'422qAS6d 




K{G{2 - 3A, l)/322tc,d/ - K(G'i(l, 1)^2(2 - A, l)))/2) 
4^) ' 



1 / 1 



'422q^CM 




K(G(2 - 3A, l)/322cacd,) 



1 1 



(87r)4 £ 



-3+-) , 
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"3 



'422qAaBd 




K(G(2-3A,l)W..) = ^i(5-^), 



'422qA6Bd 




K(G'(2 - 3A, l)/322caecd - KCdCl, 1)^2(2 - A, l))72) 



1/1 1/5 

— + 



(87r)n4£2 £ 



(4 12)) ' 



'422qA(iB6 




= K(G(2 - 3A, l)/322tcd/, + K(G(1, 1)^(1 - A, 1))J2) 



U22qAbCd 




K(G(2 - 3A, l)/322ca.cd) = " ^ + t) ' 



(J.45) 
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U22qAaac 




' 422qAabd 




'422qAbac 




' 422ciAbad 




U22ciAbbd 



l422ciAbcd — 




'422qabac 




'422qa6cd 




K(G(1, l)Gi(l, 1)^2(2 - A, 1)G(2 - 2A, 1 - A) 

-K(Gi(l, 1)6-2(2- A, l))/2) 
1 1 

K 1)2^2(2 - A, 2 - A)G(2 - 3A, 1) 



.2 
1 1 



^ = K(-G(2 - 3A, l)h22taecg) = 7^4^(3 " y 



/322w,G(2 - 3A, 1) - - X + y) 



= K(-G'(2 - 3A, l)h22tcfcg - K(G'i(l, 1)6-2(2 - A, l))/2) 

1 ^ 1 ^ 1 ^13 TT^^ 



irrY V4e2 e V 4 3 



^ = K(-G'(2 - 3A, l)/322ta,ce) = 



^ = K(G(2 - 3A, l)/322tcecc,) 



K(-G'(2 - 3A, l)l322tbgce) 



1 1/13 -2 



7r 



(87r)4£V4 3 



1 1 



(87r)4£ 



(J.46) 

The above final four-loop integrals have also been checked by using the Gegenbauer 
polynomial x-space technique (GPXT) as formulated in [H] and extended in j37lli2] . 
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J. 3. 5 Four-loop integrals with central ladder structure 



The four-loop integrals with central ladder structure and four pairwise contracted mo- 
menta in their numerators become 



M22to 




-^422t& — 



-^422tc — 




'422td 




'422te 




'422t/ 



'422tg 





K(G(2 - 3A, l)h22cadde + ^{0,(1, 1)G{2 - A, l))/2) 



4£2 £^2 8 // ' 



K(G(2-3A,l)We)=(3^J(^-^), 



K{G{2 - 3A, l)/322c6rfc. - K(G'i(l, 1)G'(2 - A, l))/2) 
1 



1 1/1 TT^ 



(87r)n 4^2 £V2 



K(-G'(2 - 3A, l)/322cadce + K(G'(l, l)G'2(2 - A, 1))72) 

7r)n 2£2 + £l^ 4>'>'' 



(Btt) 



G(2 - 3A, l)/322ccde/ = (1 - j ' 



G(2 - 3A, l)/322c6ecd - (i - y ) ' 



(J.47) 



= K(-G(2 - 3A, l)/322caecd - K(G'i(l, 1)0(1 - A, l))/2) 

(87r)n 4^2 £^4 12// 

The different ways to reexpress the four-loop integral in terms of a three-loop integral 
are based on choosing different (IR safe) external points. Cutting then out the scalar 
propagator which directly connects these points, we obtain the relations between the 
finite parts of the integrals with a central cubic vertex and with a central triangle. They 
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read 

h22cadde = h22tdffg i -^322cc/de — h22tagdf i Iz22cbdce — h22tdefg i 

h22cadce = h22tcdfg , h22ccdef — h22tafdg , -^322c6ecd = -^322tdge/ , (J -48) 

h22caecd — h22tcfdg i 

and hold only for the leading (finite) parts of the respective integrals. 

With convenient choices for the external momentum, thereby avoiding IR diver- 
gences, the required integrals with central ladder structure and six pairwise contracted 
momenta in their numerators read 




= \[K{G{1, lfG{l - A, If - 2 K(G'(1, 1)^(1 - A, l))h) 



-~ 2/42bb3a(i + h22ciabcd + 2/422ta 




(J.49) 
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J. 3. 6 Four-loop integrals for the next-to-nearest-neighbour interactions be- 
tween three internal legs 



Using the effective Feynman rules (lA.7p . the four- loop integrals which are required in 
(14.71) are determined as 



'4gq 



'4gl 




2(-l 



1 

6 



(J.50) 



4{L 



A22ciACbd — -'422qAf)Cd 



'422211 



+ 2L 



422212 — J 422213 J 



It is interesting, that /4gi is finite, as is the corresponding four-loop integral involving 
Jags- This means that integrals which involve boxes with gluon propagators do not 
contribute to the pole parts to four-loop order. It would be interesting to see if this 
continues to hold to higher loop orders. 
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